\section{仿真}
式\eqref{eqans1}和式\eqref{eqans2}都包含$\ddot{x}$和$\ddot{\theta}$，无法直接使用数值求解，否则会出现代数环。也就是说，想要计算$\ddot{x}$必须先知道$\ddot{\theta}$，而想要计算$\ddot{\theta}$也必须先知道$\ddot{x}$，形成矛盾。将式\eqref{eqans1}代入式\eqref{eqans2}，整理得到
\begin{equation}\label{eqans3}
    F=(M+m\sin^2\theta)\ddot{x}+ml\dot{\theta}^2\sin\theta-mg\sin\theta\cos\theta
\end{equation}
将式\eqref{eqans3}写成便于使用simucpp仿真的多个短方程的形式如下
\[n_1=\frac{1}{M+m\sin^2\theta}\]
\[n_2=-ml\dot{\theta}^2\sin\theta\]
\[n_3=4.9m\sin 2\theta\]
\[n_4=n_1(F+n_2+n_3)\]
\[\ddot{x}=n_4\]
式\eqref{eqans1}重新整理得到
\[n_5=\frac{1}{l}(n_4\cos\theta+9.8\sin\theta)\]
\[\ddot{\theta}=n_5\]
